Velocity correlations in dense granular flows 
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Velocity fluctuations of grains flowing down a rough inclined plane are experimentally studied. 
The grains at the free surface exhibit fluctuating motions, which are correlated over few grains 
diameters. The characteristic correlation length is shown to depend on the inclination of the plane 
and not on the thickness of the flowing layer. This result strongly supports the idea that dense 
granular flows are controlled by a characteristic length larger than the particle diameter. 

PACS numbers: 45.70.-n, 45.70.Mg, 83.50.Ax 



Dry cohesionless granular material can flow like a liq- 
uid, for example in an hourglass or on an inclined plane. 
However, the flow properties are not yet well understood 
and constitutive equations appropriate to describe the 
dense flow regime when particles are in close contact are 
still matter of debate One of the difficulty is that, 

in the dense regime, particles do not interact through bi- 
nary collisions like in dilute and agitated granular gases 
Pi, yj . They rather experience contact with several neigh- 
boring particles at the same time. The description of the 
multi-particles interactions is important in the develop- 
ment of constitutive equations. Recent theoretical mod- 
els try to take into account the existence of this contact 
network through different approaches. Some models de- 
scribe the flow as a mixture of solid and liquid phases 
or try to account for the frictional contacts by adding a 
frictional term in collisional theories 0, • Other mod- 
els consider the presence of arches in the flow 0, or 
account for some non local propagation of momentum 
through the contact network [lQ, LLl| • In this paper we 
will more specifically refer to a recent model proposed by 
Ertas and Halsey 12] to describe granular flows on an 
inclined plane. By introducing the idea of correlated mo- 
tion of grains in clusters-like structures, they have been 
able by simple scaling arguments to recover some impor- 
tant results observed for flows down inclined plane. The 
experimental study presented in this paper is inspired by 
their approach, trying to evidence the existence of corre- 
lated motions in granular flows down inclined planes. 

The inclined plane configuration is obtained when a 
granular layer of thickness h flows down a rough surface 
inclined at an angle 6 (Fig. CJ. Experiments 0, an d 
numerical simulations [lj, LUj have revealed two impor- 
tant results. First, for a given inclination a minimum 
thickness h stO p(0) exists below which no flow is possible. 
This thickness is evidenced by the deposit remaining on 
the plane once the flow stops. Second, the variation of 
the depth averaged velocity u with the inclination and 
thickness h appears to be linked to the deposit function 
hstop(6) through a flow rule: 
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where g is gravity and (3 is a constant equal to 0.13 for 
spheres. The fact that the influence of on the flow veloc- 
ity comes into play through the deposit function h stop (0) 
is not straightforward and could be simply the result of 
a coincidence 7]. However, an interesting interpretation 
of relation^has been proposed by Ertas and Halsey [l^ . 
Their idea is that motion of grains in dense granular flows 
occurs through clusters, which size is controlled by the 
stress distribution. In this picture, it is clear that no flow 
would be possible for layers thinner than the character- 
istic size of the clusters. This immediately gives an in- 
terpretation of the deposits: the deposit thickness h st0 p 
would be an indirect measurement of the cluster size. 
Based on the existence of a characteristic cluster size, 
Ertas and Halsey 12] have developed scaling arguments 
similar to a Prandtl mixing length turbulent theory and 
are able to derive the flowing ruleQ] 1]. 

This heuristic approach is appealing because it gives 
simple interpretations for the deposit thickness and the 
flow rule observed in inclined plane, but also because it 
gives new insights for the granular flow rheology. As dis- 
cussed in 1], this picture of correlated grains motion gives 
a way to link the different flow regimes observed with 
granular materials: the slow quasi-static regimes would 
correspond to strongly correlated motions whereas the 
kinetic collisional regime would correspond to a regime 
where spatial correlations vanish. To our knowledge, no 
experimental evidence exists to support this description. 

Whereas correlations are studied in details in other dis- 
ordered systems like glasses 0, or foams [HI [HI ^ , 
only few studies address the question of correlated mo- 
tions in granular flows. Most of them concern the slow 
quasi-static regime where long range correlations have 
been observed but no characteristic length has been ev- 
idenced 0, |23|. In the dense flow regime, Bonamy et 
al [2^ have investigated correlated motions in a rotat- 
ing drum. Their analysis shows that clusters of different 
sizes exist but that no characteristic length scale emerges. 
The purpose of this paper is to experimentally study spa- 
tial correlations in the inclined plane configuration and 
check if the deposit thickness observed can indeed be in- 
terpreted as an indirect measure of the correlation length 
as suggested in [l2j . 



FIG. 1: Flow threshold: (•) experimental measurements of 
the deposit h stO p(0). Inset: sketch of the experimental set up 
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The experimental set up is presented in the inset of 
Fig. [I] and is the same as the one used in A two me- 
ters long and 70 cm wide plane is made rough by gluing 
particles on it and can be inclined from horizontal at an 
angle 0. The material (glass beads 0.5 mm in diameter) 
flows from a hopper, the flow rate being controlled by the 
opening of the gate. The inclination and the thickness 
h of the layer are the two control parameters. Steady uni- 
form flows are obtained with this device in a finite range 
of inclination and thickness as shown in Fig. Q The 
flow threshold delimited by the curve h 3 top(0), has been 
determined by measuring the deposit thickness once the 
flow stops. In order to study correlated motions, we have 
chosen to precisely study the individual grains motion at 
the free surface using a high speed video camera placed 
above the plane. From movies recorded between 500 and 
1000 frames per second, we measure the velocity of the 
particles using a precise particle tracking algorithm. The 
mean free surface velocity U s e x is computed by averaging 
over 1 second of flow and over all the detected particles. 
It is then subtracted to all individual particle velocities to 
get the instantaneous fluctuating velocity field. For each 
time t of our movie, we then have the position (x l ,y 1 } 
of particle z, and its velocity relative to the mean free 
surface motion (u x ,u y ). Example of fluctuating velocity 
fields are presented in Figs. [21 (a) and (b) for two differ- 
ent inclinations. Fig. [3 (a) obtained at low inclination 
clearly exhibits correlated motions. By contrast, the ve- 
locity field is much more disordered with less correlations 
at high inclination as shown in Fig. 0(b). Before study- 
ing spatial correlations, we first analyze the amplitude 
of the fluctuations by computing the mean fluctuating 
velocity 5V = (< (u x ) 2 + (u y ) 2 >) 1/2 averaged over the C xx (Ax) = 
particles and over lsecond of flow. In Fig. the dimen- 
sionless fluctuating velocity amplitude SVj yfgd is plotted 
as a function of the mean shear rate i.e. U s /h^/ g / d. All 
the data obtained for different inclinations and different 
flow thicknesses collapse on a single curve, showing that 
the fluctuations at the free surface are entirely driven by 
the shear rate. The plot seems to indicate an affine re- 
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FIG. 2: Fluctuating velocity of particles at the free surface 
(a) = 21°, h/d = 13.7;(b) = 26°, h/d=9.3 . 



lation between fluctuations and shear rate. However, if 
one tries to fit the data by a power law in order to com- 
pare with results published in the literature, one get a 
power equal to 0.7 , compatible with measurements in 
other flow configurations 

We now turn to the study of the spatial correlations. 
To this end, we defined four correlations function C a p 
corresponding to the correlation of the (3 velocity compo- 
nent along the a direction. For example, the correlation 
C xx (Ax) at a distance (Ax) is defined as follow: 

Et (Eij + Ax - - y*) 

Et (Eij S(x i + Ax - xJ)5(y* - yi)) 

where the function S(x) is a peaked function. We have 
chosen a Gaussian of width 0.4 particle diameter for S 
but have checked that the results are insensitive to this 
choice, as long as 5 is narrow. Fig. 21 shows the four 
correlation functions for = 20° and h/d = 19.8. The 
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four functions have the same characteristics: a first rapid 
decrease, followed after 2d by a roughly exponential de- 
crease over a characteristic length scale of few particles 
diameters. The function C xx is larger than the others, 
meaning that the correlations are stronger for the longi- 
tudinal velocity along the flow direction. In the following, 
we show the influence of the inclination and flow thick- 
ness h on C xx only, but the same qualitative behaviors 
hold for the three other correlation functions. 




FIG. 4: The four normalized correlation functions for — 20° 
and h/d = 19.8; C xx (A) , C yy (O) , C yx (▼), C xy ( □). 

The three plots in Fig. present the normalized cor- 
relation function C xx /C xx (0) at three different inclina- 
tions for different flow thicknesses h. The striking result 
shown by these figures is the independence of the corre- 
lation function with the thickness h. At a given angle 
0, the normalized correlation function remains identical 
when one increases the flow thickness. This means that 
the spatial range of correlation is not controlled by the 
thickness but varies with inclination only. As a conse- 



quence, to each inclination corresponds a single correla- 
tion function, which can be computed by averaging the 
correlation functions obtained for different thicknesses. 
The results obtained for all the inclinations are summa- 
rized in Fig. El The correlation function is steeper and 
steeper when increasing the plane inclination 0, show- 
ing that the correlation length decreases. In order to 
quantify this effect, we define a correlation length L xx as 
the length where the correlation is equal to 0.07 times 
its value at the origin: C XX (L XX ) = 0.07C XX (0). Other 
definition gives similar qualitative results. Our measure- 
ments then show that L xx is a decreasing function of 
the inclination as shown in Fig. At high inclina- 
tion, close to the maximum angle where steady uniform 
flows are observed, the correlation length is close to one 
particle diameter, which means that almost no correla- 
tion between neighboring particles exists. However, the 
correlation length dramatically increases for low inclina- 
tions when approaching the minimum angle at which a 
flow is possible. The same evolution with inclination, 
although less pronounced, is observed for the four corre- 
lation lengths L a p (Fig. 0. This variation of the correla- 
tion length with inclination is reminiscent of the deposit 
function h stO p(0) (Fig. [lj. The deposit seems to diverge 
at the minimum angle where flow is possible and goes to 
zero at higher inclination. Both quantities L xx and h st0 p 
thus varies in the same way with inclination as shown 
in inset of Fig. [3 where L xx /d is plotted as a function 
hstop/d for the different inclinations. 

Given the error bars in the measurement of the corre- 
lation length, it is difficult to conclude about a precise 
relation between L xx and h s t op - However, despite the 
uncertainties, three major conclusions can be drawn: 
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FIG. 6: Mean normalized correlation function for different 0. 



when varying the roughness. An interesting question is 
how the correlation length would be also affected. De- 
spite these open questions, this work gives a first exper- 
imental evidence to the proposition by Ertas and Halsey 
[l^ that the flow threshold h s t op (0) observed in the in- 
clined plane configuration is the signature of correlated 
motions. More generally speaking, our results strongly 
support the idea that the rheology of granular material 
in the dense regime might be controlled by correlated mo- 
tion of grains on a length scale larger than the particle 
diameter. The microscopic origin of the spatial corre- 
lations and how to properly take them into account in 
writing constitutive laws, remain important challenges. 

This work has benefited from fruitful discussions with 
B. Andreotti, O. Dauchot, Y. Forterre and P. Jop. 
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FIG. 7: The four correlation lengths L xx (a) , L yy (O) , L yx 
(T), L xy ( □)., as function of 0; inset: L xx as a function of 
h s top- The error bars give the measurement dispersion when 
the correlation lengths are computed from individual runs and 
not from the mean correlation functions of Fig. |S] 



i) Correlated motions exist at the free surface of a dense 
granular flow on a rough inclined plane. 

ii) The spatial extend of the correlated motion is not 
controlled by the thickness of the flowing layer but is 
fixed by the inclination only. 

iii) The correlation length is maximum at low incli- 
nation and decreases at high inclination, in a similar 
way the minimum thickness for which flow is possible 
decreases with inclination. 

These conclusions only hold for the particles at the free 
surface. It would be interesting to know to which extend 
they still hold in the bulk of the flow, a study that could 
be possible in numerical simulations. The study also 
only concerns the generic case where the roughness of the 
plane is made of the same beads as in the bulk. However, 
a recent study Q shows that the flow properties varies 
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